Teleost-specific Whole Genome Duplication Research Articles

Comparative genomics studies on the stk gene family in vertebrates: From the bighead carp (Hypophthalmichthys nobilis) genome

The mammalian sterile 20-like (MST) family belongs to the serine/threonine protein kinase (STK) superfamily and participates in a variety of biological processes, such as cell apoptosis, polarity, migration, immune regulation, inflammatory responses, and cancer. In the economically important bighead carp (Hypophthalmichthys nobilis), the STK gene family and immune-related biological functions may be helpful in increasing its economic yield. However, the comprehensive role of STKs in the bighead carp remains unclear. In this study, the five stk sequences from the bighead carp were divided into two classes: stk3/4 and stk24/25/26. Gene structure and motif prediction analyses confirmed that stk is conserved in the bighead carp. Compared to 26 other vertebrate species, teleosts (including bighead carp) possess more stk members because of teleost-specific whole-genome duplication. Synteny analysis revealed that stk3, stk24, stk25, and stk26 have been relatively conserved in bighead carp during evolution. Meanwhile, stk4 was lost in most Cyprinid species, including bighead carp, during evolution. RNA-seq data revealed that STK expression was associated with various pathogens, and the expression of these STKs (Hnstk3, Hnstk24a, Hnstk24b, Hnstk25, and Hnstk26) was different in seven tissues of bighead carp. In addition, we showed that STK expression levels were dramatically altered in the head kidney and that stk24 was involved in defense against Aeromonas hydrophila. This study provides a molecular basis for the analysis of stk function in bighead carp, and can be used as a reference for further phylogenomics.

Novel Dimeric Architecture of an IFN-γ-Related Cytokine Provides Insights into Subfunctionalization of Type II IFNs in Teleost Fish.

Gene duplication leads to subfunctionalization of paralogs. In mammals, IFN-γ is the sole member of the type II IFN family and binds to a receptor complex consisting of IFN-γR1 and IFN-γR2. In teleost fish, IFN-γ and its receptors have been duplicated due to the teleost-specific whole-genome duplication event. In this study, the functions of an IFN-γ-related (IFN-γrel) cytokine were found to be partially retained relative to IFN-γ in grass carp (Ctenopharyngodon idella [CiIFN-γrel]). CiIFN-γrel upregulated the expression of proinflammatory genes but had lost the ability to activate genes involved in Th1 response. The results suggest that CiIFN-γrel could have been subfunctionalized from CiIFN-γ. Moreover, CiIFN-γrel induced STAT1 phosphorylation via interaction with duplicated homologs of IFN-γR1 (cytokine receptor family B [CRFB] 17 and CRFB13). Strikingly, CiIFN-γrel did not bind to the IFN-γR2 homolog (CRFB6). To gain insight into the subfunctionalization, the crystal structure of CiIFN-γrel was solved at 2.26 Å, revealing that it forms a homodimer that is connected by two pairs of disulfide bonds. Due to the spatial positions of helix A, loop AB, and helix B, CiIFN-γrel displays a unique topology that requires elements from two identical monomers to form a unit that is similar to IFN-γ. Further, mutagenesis analyses identified key residues interacting with CiIFN-γrel receptors and those required for the biological functions. Our study can help understand the subfunctionalization of duplicated IFN-γ paralogs in fish.

Three IRF4 paralogs act as negative regulators of type Ⅰ IFN responses in yellow catfish (Pelteobagrus fulvidraco)

IRF4 is a master member of the interferon regulatory factor (IRF) family playing vital regulatory roles in immune system development and function. Tetrapods have a single-copy IRF4 gene, while teleosts harbor duplicated IRF4 genes. This work describes three IRF4 paralogs from yellow catfish (Pelteobagrus fulvidraco), designated PfIRF4A, PfIRF4B and PfIRF4B-like. These genes all contain a typical IRF structural architecture. Phylogenic and synteny analyses indicate that they should arise from the teleost-specific whole-genome duplication. PfIRF4 genes are abundantly expressed in the immune-related tissues and upregulated by PolyI:C, LPS, and Edwardsiella ictaluri. Ectopic expression of these genes inhibits the activation of fish type Ⅰ IFN promoters and downregulates the transcription levels of IFN-responsive genes, thus allowing the efficient replication of a fish rhabdovirus, spring viremia of carp virus (SVCV). PfIRF4s possess a repressive effect on MyD88-mediated activation of IFN and NF-κB. Some differences are observed between each individual paralog. PfIRF4B is the main form expressed across the tissues and the most up-regulated one after pathogen induction. It exerts a stronger inhibitory effect on IFN antiviral response than the other two paralogs. PfIRF4A and PfIRF4B-like are primarily present in the nucleus, while PfIRF4B displays colocalization and direct associations with MyD88 in the cytoplasm. Overall, the data demonstrate that three PfIRF4 paralogs show shared and individual functional properties in the negative regulation of type Ⅰ IFN response.

Auts2 Features and Expression Are Highly Conserved during Evolution Despite Different Evolutionary Fates Following Whole Genome Duplication.

The AUTS2 gene plays major roles during brain development and is associated with various neuropathologies including autism. Data in non-mammalian species are scarce, and the aim of our study was to provide a comprehensive analysis of auts2 evolution in teleost fish, which are widely used for in vivo functional analysis and biomedical purposes. Comparative genomics in 78 species showed that auts2a and auts2b originate from the teleost-specific whole genome duplication (TGD). auts2a, which is highly similar to human AUTS2, was almost systematically retained following TGD. In contrast, auts2b, which encodes for a shorter protein similar to a short human AUTS2 isoform, was lost more frequently and independently during evolution. RNA-seq analysis in 10 species revealed a highly conserved profile with predominant expression of both genes in the embryo, brain, and gonads. Based on protein length, conserved domains, and expression profiles, we speculate that the long human isoform functions were retained by auts2a, while the short isoform functions were retained by auts2a and/or auts2b, depending on the lineage/species. auts2a showed a burst in expression during medaka brain formation, where it was expressed in areas of the brain associated with neurodevelopmental disorders. Together, our data suggest a strong conservation of auts2 functions in vertebrates despite different evolutionary scenarios in teleosts.

New Insights Into the Evolution of Corticotropin-Releasing Hormone Family With a Special Focus on Teleosts.

Corticotropin-releasing hormone (CRH) was discovered for its role as a brain neurohormone controlling the corticotropic axis in vertebrates. An additional crh gene, crh2, paralog of crh (crh1), and likely resulting from the second round (2R) of vertebrate whole genome duplication (WGD), was identified in a holocephalan chondrichthyan, in basal mammals, various sauropsids and a non-teleost actinopterygian holostean. It was suggested that crh2 has been recurrently lost in some vertebrate groups including teleosts. We further investigated the fate of crh1 and crh2 in vertebrates with a special focus on teleosts. Phylogenetic and synteny analyses showed the presence of duplicated crh1 paralogs, crh1a and crh1b, in most teleosts, resulting from the teleost-specific WGD (3R). Crh1b is conserved in all teleosts studied, while crh1a has been lost independently in some species. Additional crh1 paralogs are present in carps and salmonids, resulting from specific WGD in these lineages. We identified crh2 gene in additional vertebrate groups such as chondrichthyan elasmobranchs, sarcopterygians including dipnoans and amphibians, and basal actinoperygians, Polypteridae and Chondrostei. We also revealed the presence of crh2 in teleosts, including elopomorphs, osteoglossomorphs, clupeiforms, and ostariophysians, while it would have been lost in Euteleostei along with some other groups. To get some insights on the functional evolution of the crh paralogs, we compared their primary and 3D structure, and by qPCR their tissue distribution, in two representative species, the European eel, which possesses three crh paralogs (crh1a, crh1b, crh2), and the Atlantic salmon, which possesses four crh paralogs of the crh1-type. All peptides conserved the structural characteristics of human CRH. Eel crh1b and both salmon crh1b genes were mainly expressed in the brain, supporting the major role of crh1b paralogs in controlling the corticotropic axis in teleosts. In contrast, crh1a paralogs were mainly expressed in peripheral tissues such as muscle and heart, in eel and salmon, reflecting a striking subfunctionalization between crh1a and b paralogs. Eel crh2 was weakly expressed in the brain and peripheral tissues. These results revisit the repertoire of crh in teleosts and highlight functional divergences that may have contributed to the differential conservation of various crh paralogs in teleosts.

CRISPR Knockouts of pmela and pmelb Engineered a Golden Tilapia by Regulating Relative Pigment Cell Abundance.

Premelanosome protein (pmel) is a key gene for melanogenesis. Mutations in this gene are responsible for white plumage in chicken, but its role in pigmentation of fish remains to be demonstrated. In this study, we found that most fishes have 2 pmel genes arising from the teleost-specific whole-genome duplication. Both pmela and pmelb were expressed at high levels in the eyes and skin of Nile tilapia. We mutated both genes in tilapia using CRISPR/Cas9. Homozygous mutation of pmela resulted in yellowish body color with weak vertical bars and a hypopigmented retinal pigment epithelium (RPE) due to significantly reduced number and size of melanophores. In contrast, we observed an increased number and size of xanthophores in mutants compared to wild-type fish. Homozygous mutation of pmelb resulted in a similar, but milder phenotype than pmela-/- mutants. Double mutation of pmela and pmelb resulted in loss of additional melanophores compared to the pmela-/- mutants, and also an increase in the number and size of xanthophores, producing a golden body color. The RPE pigmentation of pmela-/-;pmelb-/- was similar to pmela-/- mutants, with much less pigmentation than pmelb-/- mutants and wild-type fish. Taken together, our results indicate that, although both pmel genes are important for the formation of body color in tilapia, pmela plays a more important role than pmelb. To our knowledge, this is the first report on mutation of pmelb or both pmela;pmelb in fish. Studies on these mutants suggest new strategies for breeding golden tilapia, and also provide a new model for studies of pmel function in vertebrates.

Evolution of the Parvalbumin Genes in Teleost Fishes after the Whole-Genome Duplication

Parvalbumin is considered a major fish allergen. Here, we report the molecular evolution of the parvalbumin genes in bony fishes based on 19 whole genomes and 70 transcriptomes. We found unexpectedly high parvalbumin diversity in teleosts; three main gene types (pvalb-α, pvalb-β1, and pvalb-β2, including oncomodulins) originated at the onset of vertebrates. Teleosts have further multiplied the parvalbumin gene repertoire up to nine ancestral copies—two copies of pvalb-α, two copies of pvalb-β1, and five copies of pvalb-β2. This gene diversity is a result of teleost-specific whole-genome duplication. Two conserved parvalbumin genomic clusters carry pvalb-β1 and β2 copies, whereas pvalb-α genes are located separately in different linkage groups. Further, we investigated parvalbumin gene expression in 17 tissues of the common carp (Cyprinus carpio), a species with 21 parvalbumin genes in its genome. Two pvalb-α and eight pvalb-β2 copies are highly expressed in the muscle, while two alternative pvalb-α copies show expression in the brain and the testes, and pvalb-β1 is dominant in the retina and the kidney. The recent pairs of muscular pvalb-β2 genes show differential expression in this species. We provide robust genomic evidence of the complex evolution of the parvalbumin genes in fishes.

Evolution of inwardly rectifying potassium channels and their gene expression in zebrafish embryos.

Inwardly rectifying potassium channels are essential for normal potassium homeostasis, maintaining the cellular resting membrane potential, and regulating electrolyte transportation. Mutations in Kir channels have been known to cause debilitating diseases ranging from neurological abnormalities to renal and cardiac failures. Many efforts have been made to understand their protein structures, physiological functions, and pharmacological modifiers. However, their expression and functions during embryonic development remain largely unknown. Using zebrafish as a model, we identified and renamed 31 kir genes. We also analyzed Kir gene evolution by phylogenetic and syntenic analyses. Our data indicated that the four subtypes of the Kir genes might have already evolved out in chordates. These vertebrate Kir genes most likely resulted from both whole-genome duplications and tandem duplications. In addition, we examined zebrafish kir gene expression during early embryogenesis. Each subgroup's genes showed similar but distinct gene expression domains. The gene expression of ohnologous genes from teleost-specific whole-genome duplication indicated subfunctionalization. Varied temporal gene expression domains suggest that Kir channels may be needed for embryonic patterning or regulation. Our phylogenetic and developmental analyses of Kir channels shed light on their evolutionary history and potential functions during embryogenesis related to congenital diseases and human channelopathies.

Partitioning of gene expression among zebrafish photoreceptor subtypes

Vertebrate photoreceptors are categorized into two broad classes, rods and cones, responsible for dim- and bright-light vision, respectively. While many molecular features that distinguish rods and cones are known, gene expression differences among cone subtypes remain poorly understood. Teleost fishes are renowned for the diversity of their photoreceptor systems. Here, we used single-cell RNA-seq to profile adult photoreceptors in zebrafish, a teleost. We found that in addition to the four canonical zebrafish cone types, there exist subpopulations of green and red cones (previously shown to be located in the ventral retina) that express red-shifted opsin paralogs (opn1mw4 or opn1lw1) as well as a unique combination of cone phototransduction genes. Furthermore, the expression of many paralogous phototransduction genes is partitioned among cone subtypes, analogous to the partitioning of the phototransduction paralogs between rods and cones seen across vertebrates. The partitioned cone-gene pairs arose via the teleost-specific whole-genome duplication or later clade-specific gene duplications. We also discovered that cone subtypes express distinct transcriptional regulators, including many factors not previously implicated in photoreceptor development or differentiation. Overall, our work suggests that partitioning of paralogous gene expression via the action of differentially expressed transcriptional regulators enables diversification of cone subtypes in teleosts.

Gene clusters related to metamorphosis in Solea senegalensis are highly conserved

The flatfish, Solea senegalensis has considerable scientific interest and commercial value. The metamorphosis in this species occurs between 12 and 19 days after hatching and it takes about 1 week to complete. Eleven Bacterial Artificial Chromosomes (BAC) clones containing the various candidate genes involved in the process of metamorphosis: thyroxine 5 deiodinase 3 (dio3); forkhead box protein E4 (foxe4); melatonin receptor type 1C (mel1c); calsequestrin 1b (casq1b); thyrotropin subunit beta (tshβ); thyrotropin-releasing hormone receptor 1, 2, and 3 (trhr1, trhr2, trhr3); thyroid hormone receptor α a and b (thrαa, thrαb); and thyroid hormone receptor beta (thrβ) were analyzed by multiple Fluorescence in situ Hybridization (mFISH) and Next Generation Sequencing (NGS) techniques. The mFISH technique localized the 11 BAC clones on 12 different chromosome pairs because three of them, specifically the trhr1a, trhr2 and thrβ BAC clones, showed double signals. This signal duplication indicates a duplication of the genomic region inserted within the BAC clone, which provides evidence for the Teleost-Specific Whole Genome Duplication (TS-WGD). Micro-synteny and phylogenetic analysis showed that Cynoglossus semilaevis is the nearest species to S. senegalensis and that Danio rerio is the most distant one. The tshβ BAC clone was highly conserved as the genes belonging to this BAC were located on a single chromosome in all the species studied. These genes participate in proliferation, migration and cell-death, which are key processes during metamorphosis. Overall, micro-synteny analysis showed that most candidate genes are found in conserved genomic surroundings.

Visual pigment evolution in Characiformes: The dynamic interplay of teleost whole‐genome duplication, surviving opsins and spectral tuning

Vision represents an excellent model for studying adaptation, given the genotype-to-phenotype map that has been characterized in a number of taxa. Fish possess a diverse range of visual sensitivities and adaptations to underwater light, making them an excellent group to study visual system evolution. In particular, some speciose but understudied lineages can provide a unique opportunity to better understand aspects of visual system evolution such as opsin gene duplication and neofunctionalization. In this study, we showcase the visual system evolution of neotropical Characiformes and the spectral tuning mechanisms they exhibit to modulate their visual sensitivities. Such mechanisms include gene duplications and losses, gene conversion, opsin amino acid sequence and expression variation, and A1 /A2 -chromophore shifts. The Characiforms we studied utilize three cone opsin classes (SWS2, RH2, LWS) and a rod opsin (RH1). However, the characiform's entire opsin gene repertoire is a product of dynamic evolution by opsin gene loss (SWS1, RH2) and duplication (LWS, RH1). The LWS- and RH1-duplicates originated from a teleost specific whole-genome duplication as well as characiform-specific duplication events. Both LWS-opsins exhibit gene conversion and, through substitutions in key tuning sites, one of the LWS-paralogues has acquired spectral sensitivity to green light. These sequence changes suggest reversion and parallel evolution of key tuning sites. Furthermore, characiforms' colour vision is based on the expression of both LWS-paralogues and SWS2. Finally, we found interspecific and intraspecific variation in A1 /A2 -chromophores proportions, correlating with the light environment. These multiple mechanisms may be a result of the diverse visual environments where Characiformes have evolved.

Molecular characterization and expression analysis of foxo3l in response to exogenous hormones in black rockfish (Sebastes schlegelii).

As a crucial member of the Forkhead Box family, class O (FoxO) plays an essential role in growth, cell differentiation, metabolism, immunization, and apoptosis. Meanwhile, FoxO3 is the primary regulator and effective inhibitor of primordial follicle activation. In this study, seven foxo genes were identified in black rockfish (Sebastes schlegelii), including two foxo1 genes (foxo1a, foxo1b), two foxo3 genes (foxo3, foxo3l), one foxo4 gene, and two foxo6 genes (foxo6a, foxo6b). foxo3l was derived from teleost-specific whole-genome duplication events. Evaluation of tissue expression pattern revealed that foxo3l displayed sexually dimorphic expression with a high level in the ovary and spatial expression only in the cytoplasm of follicle cells and oocytes. When the ovaries were stimulated by estrogen and gonadotropin, foxo3l expression was remarkably reduced, and the effect of androgen was completely different. We considered that foxo3l lost its ability to inhibit follicular precocity because of mass ovulation by hormone stimulation, resulting in its decreased expression. Such evidence indicated that foxo3l is an important regulator of reproduction-related functions in black rockfish. This study provides new insights into foxo3l genes for further functional research in teleost.

Beta-Hydroxysteroid Dehydrogenase Genes in Orange-Spotted Grouper (Epinephelus coioides): Genome-Wide Identification and Expression Analysis During Sex Reversal.

Beta-hydroxysteroid dehydrogenases (β-HSDs) are a group of steroidogenic enzymes that are involved in steroid biosynthesis and metabolism, and play a crucial role in mammalian physiology and development, including sex determination and differentiation. In the present study, a genome-wide analysis identified the numbers of β-hsd genes in orange-spotted grouper (Epinephelus coioides) (19), human (Homo sapiens) (22), mouse (Mus musculus) (24), chicken (Gallus gallus) (16), xenopus (Xenopus tropicalis) (24), coelacanth (Latimeria chalumnae) (17), spotted gar (Lepisosteus oculatus) (14), zebrafish (Danio rerio) (19), fugu (Takifugu rubripes) (19), tilapia (Oreochromis niloticus) (19), medaka (Oryzias latipes) (19), stickleback (Gasterosteus aculeatus) (17) and common carp (Cyprinus carpio) (27) samples. A comparative analysis revealed that the number of β-hsd genes in teleost fish was no greater than in tetrapods due to gene loss followed by a teleost-specific whole-genome duplication event. Based on transcriptome data from grouper brain and gonad samples during sex reversal, six β-hsd genes had relatively high expression levels in the brain, indicating that these genes may be required for neurogenesis or the maintenance of specific biological processes in the brain. In the gonad, two and eight β-hsd genes were up- and downregulated, respectively, indicating their important roles in sex reversal. Our results demonstrated that β-hsd genes may be involved in the sex reversal of grouper by regulating the synthesis and metabolism of sex steroid hormones.

Identification and Characterization of Carnitine Palmitoyltransferase 1 (cpt 1) Genes in Nile Tilapia, Oreochromis niloticus.

Four cpt 1 genes (cpt 1α1a, cpt 1α2a, cpt 1α2b, and cpt 1β) were identified in the Nile tilapia genome. Two transmembrane helix domains (TMH) were identified for Cpt 1α1a, Cpt 1α2a, and Cpt 1β, while Cpt 1α2b had only one TMH domain. Evidence was found of conserved gene synteny between cpt 1 genes from Nile tilapia and the cpt 1/CPT 1 genes of zebrafish and human. Phylogenetic analysis showed that Nile tilapia Cpt 1 sequences clustered in distinct clades with their orthologous Cpt 1/CPT 1 from other vertebrates. Nile tilapia cpt 1α1a, cpt 1α2a, cpt 1α2b, and cpt 1β contain 18 coding exons encoding polypeptides of 771, 784, 788, and 786 amino acids in length, respectively. The cpt 1 genes were determined in all the tested tissues with varying tissue distribution patterns. These findings suggest that (1) cpt 1α1a, cpt 1α2a, and cpt 1α2b arose in the Nile tilapia genome as a result of the teleost-specific whole-genome duplication; (2) nonfunctionalization is the most likely cause of the loss of cpt 1α1b in the Nile tilapia genome; (3) the different tissue-specific transcription of cpt 1α2a and cpt 1α2b may be either due to the sub- or the neo-functionalization of transcriptional control side.

Genome-Wide Identification and Expression Profile of the Sox Gene Family During Embryo Development in Large Yellow Croaker, Larimichthys crocea.

Members of the Sox gene family play crucial roles during reproduction and development, but their genome-wide identification has not yet been performed in large yellow croaker, Larimichthys crocea. In this study, a total of 26 members of the Sox gene family were identified from the genome of large yellow croaker and classified into seven subgroups based on the conserved HMG-box domain they contain. Among the identified Sox gene family members, eight belonged to the SoxB subgroup (five in B1 and three in B2), four belonged to the SoxC subgroup, four belonged to the SoxD subgroup, six belonged to the SoxE subgroup, three belonged to the SoxF subgroup, and one belonged to the SoxK subgroup. During evolution, members of the SoxE subgroup (Sox8, Sox9, Sox10), Sox1, Sox4, Sox6, and Sox11 evolved into two copies, which may be a result of teleost-specific whole-genome duplication. Sox genes were distributed unevenly across 15 chromosomes. The number of introns in large yellow croaker Sox genes varied from 0 to 14. Results of the expression profile during embryogenesis revealed that most of the members of the Sox gene family had lower expression, except several Sox genes, and expression patterns also differed among each Sox gene group and duplicated gene. This study systematically characterized and analyzed the Sox gene family in large yellow croaker and provided new insights into its function during embryogenesis.

Evolutionary significance and regulated expression of Tdrd family genes in gynogenetic Japanese flounder (Paralichthys olivaceus)

Tudor domain-containing proteins (TDRDs) are highly conserved among organisms and have a function in gonads to regulate gametogenesis and genome stability through the piwi-interacting RNA (piRNA) pathway. With diverse sexual development patterns in teleost species, the evolution and function of Tdrd genes among teleosts remain unclear. Here, we identified and characterized 12 Tdrd genes (PoTdrds) in Japanese flounder (Paralichthys olivaceus) which represents dramatic sexual dimorphic metrics and sex reversal during sex differentiation. Phylogenetic and comparative synteny indicated the gain and loss of Tdrd genes after teleost-specific whole-genome duplication (3R-WGD). Tdrd1, Tdrd5, Tdrd6 and Ecat8 were abundantly expressed in their gonads. Four PoTdrds (Tdrd6, Tdrd7b, Tdrd9 and Ecat8) represented significant male-biased expression in gynogenetic and wild-type Japanese flounder gonads (p < .01). This finding indicated their important roles in spermatogenesis of P. olivaceus. Some PoTdrds were either highly up-regulated in gynogenetic testis (Tdrd3, Tdrd5, Tdrd7b and Ecat8) or down-regulated in gynogenetic ovary (Tdrkh, Tdrd3, Tdrd6l) compared with wild-type gonads (p < .05). Molecular evolution tests revealed that the selective pressure of Tdrd6/6l differed between ancestral aquatic and terrestrial organisms with 13 positively selected sites found in the ancestral lineages of teleost Tdrd6. Expression profile analysis suggested that PoTdrd6 differed significantly from PoTdrd6l, indicating sub-functionalization after 3R-WGD. All these results are important for the functional annotation of Tdrd genes and can benefit the further deciphering of Tdrd functions during gonadal development and gametogenesis of teleost fish.

Characterization of paralogous uncx transcription factor encoding genes in zebrafish.

The paired-type homeodomain transcription factor Uncx is involved in multiple processes of embryogenesis in vertebrates. Reasoning that zebrafish genes uncx4.1 and uncx are orthologs of mouse Uncx, we studied their genomic environment and developmental expression. Evolutionary analyses indicate the zebrafish uncx genes as being paralogs deriving from teleost-specific whole-genome duplication. Whole-mount in situ mRNA hybridization of uncx transcripts in zebrafish embryos reveals novel expression domains, confirms those previously known, and suggests sub-functionalization of paralogs. Using genetic mutants and pharmacological inhibitors, we investigate the role of signaling pathways on the expression of zebrafish uncx genes in developing somites. In identifying putative functional role(s) of zebrafish uncx genes, we hypothesized that they encode transcription factors that coordinate growth and innervation of somitic muscles.

Oxygenation properties of hemoglobin and the evolutionary origins of isoform multiplicity in an amphibious air-breathing fish, the blue-spotted mudskipper (Boleophthalmus pectinirostris).

Among the numerous lineages of teleost fish that have independently transitioned from obligate water breathing to facultative air breathing, evolved properties of hemoglobin (Hb)-O2 transport may have been shaped by the prevalence and severity of aquatic hypoxia (which influences the extent to which fish are compelled to switch to aerial respiration) as well as the anatomical design of air-breathing structures and the cardiovascular system. Here, we examined the structure and function of Hbs in an amphibious, facultative air-breathing fish, the blue-spotted mudskipper (Boleophthalmus pectinirostris). We also characterized the genomic organization of the globin gene clusters of the species and we integrated phylogenetic and comparative genomic analyses to unravel the duplicative history of the genes that encode the subunits of structurally distinct mudskipper Hb isoforms (isoHbs). The B. pectinirostris isoHbs exhibit high intrinsic O2 affinities, similar to those of hypoxia-tolerant, water-breathing teleosts, and remarkably large Bohr effects. Genomic analysis of conserved synteny revealed that the genes that encode the α-type subunits of the two main adult isoHbs are members of paralogous gene clusters that represent products of the teleost-specific whole-genome duplication. Experiments revealed no appreciable difference in the oxygenation properties of co-expressed isoHbs in spite of extensive amino acid divergence between the alternative α-chain subunit isoforms. It therefore appears that the ability to switch between aquatic and aerial respiration does not necessarily require a division of labor between functionally distinct isoHbs with specialized oxygenation properties.

Hypoxia-inducible transcription factors in fish: expression, function and interconnection with the circadian clock.

The hypoxia-inducible transcription factors are key regulators for the physiological response to low oxygen availability. In vertebrates, typically three Hif-α isoforms, Hif-1α, Hif-2α and Hif-3α, are expressed, each of which, together with Hif-1β, may form a functional heterodimer under hypoxic conditions, controlling expression of hundreds of genes. A teleost-specific whole-genome duplication complicates the analysis of isoform-specific functions in fish, but recent studies suggest that the existence of paralogues of a specific isoform opens up the possibility for a subfunctionalization. In contrast to during development inside the uterus, fish eggs are freely accessible and studies analyzing Hif expression in fish embryos during development have revealed that Hif proteins are not only controlling the hypoxic response, but are also crucial for proper development and organ differentiation. Significant advances have been made in our knowledge about tissue-specific functions of Hif proteins, especially with respect to gill or gonadal tissue. The hypoxia signalling pathway is known to be tightly and mutually intertwined with the circadian clock in zebrafish and mammals. Recently, a mechanistic explanation for the hypoxia-induced dampening of the transcriptional clock was detected in zebrafish, including also metabolically induced alterations of cellular redox signalling. In turn, MAP kinase-mediated H2O2 signalling modulates the temporal expression of Hif-1α protein, similar to the redox regulation of the circadian clock itself. Once again, the zebrafish has emerged as an excellent model organism with which to explore these specific functional aspects of basic eukaryotic cell biology.

Centromere evolution and CpG methylation during vertebrate speciation

Centromeres and large-scale structural variants evolve and contribute to genome diversity during vertebrate speciation. Here, we perform de novo long-read genome assembly of three inbred medaka strains that are derived from geographically isolated subpopulations and undergo speciation. Using single-molecule real-time (SMRT) sequencing, we obtain three chromosome-mapped genomes of length ~734, ~678, and ~744Mbp with a resource of twenty-two centromeric regions of length 20–345kbp. Centromeres are positionally conserved among the three strains and even between four pairs of chromosomes that were duplicated by the teleost-specific whole-genome duplication 320–350 million years ago. The centromeres do not all evolve at a similar pace; rather, centromeric monomers in non-acrocentric chromosomes evolve significantly faster than those in acrocentric chromosomes. Using methylation sensitive SMRT reads, we uncover centromeres are mostly hypermethylated but have hypomethylated sub-regions that acquire unique sequence compositions independently. These findings reveal the potential of non-acrocentric centromere evolution to contribute to speciation.